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Abstract 

We prove existence of strongly continuous evolution systems in L 2 for 
Schrodinger-type equations with non-Lipschitz coefficients in the princi- 
pal part. The underlying operator structure is motivated from models 
of paraxial approximations of wave propagation in geophysics. Thus, 
the evolution direction is a spatial coordinate (depth) with additional 
pseudodifferential terms in time and low regularity in the lateral space 
variables. We formulate and analyze the Cauchy problem in distribu- 
tion spaces with mixed regularity. The key point in the evolution system 
construction is an elliptic regularity result, which enables us to precisely 
determine the common domain of the generators. The construction of a 
solution with low regularity in the coefficients is the basis for an inverse 
analysis which allows to infer the lack of lateral regularity in the medium 
from measured data. 
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1 Introduction 



The paraxial equations in models of wave propagation are based on parabolic 
symbol approximations in theories of wave operators. They have been exten- 
sively applied in integrated optics, underwater acoustic tomography as well as 
reflection seismic imaging (cf. [21 [H]). They have also entered the analysis of 
time-reversal mirror experiments with waves taking into account stochastic vari- 
ations in the wave speed (cf. [2]). The paraxial equations can be derived from the 
reduced wave or Helmholtz equation and, since they split the wave fields accord- 
ing to a prescribed principal direction of propagation, are also called one-way 
wave equations. In particular, the leading-order parabolic symbol approxima- 
tion, also referred to as the narrow-angle or beam-propagation approximation, 
leads to model equations of Schrodinger-type. The well-posedness of the one- 
way wave Cauchy problems has been discussed by [6l [22] . The methodologies 
developed to date, however, have assumed smoothness of the wave speed func- 
tion (i.e., the coefficients in the wave operators). 

In the analysis presented here, we depart from this smoothness assumption by 
allowing the coefficients to be of any Holder regularity between zero and one, 
but typically non-Lipschitz. The existence of distributional solutions to second 
order strictly hyperbolic equations in general may fail below Log-Lipschitz reg- 
ularity of the coefficients (cf. 4!). In case of Holder regularity 2 or higher a 
constructive approach for hyperbolic evolution equations has been developed in 
PP, thereby extending results on propagation of singularities. The particular 
equation considered here is derived from such a second order equation, but the 
existence of its solution is not restricted by the same conditions. Indeed, we 
exploit the framework of Sobolev space techniques, in particular, multiplication 
of distributions in scales of Sobolev spaces, to construct a strongly continuous 
evolution system. The novelty of the paper lies in the method of construction 
which not only provides a solution concept for the paraxial wave equation with 
low coefficient regularity, but also allows us to investigate how the coefficient 
regularity influences the solution. 

The class of coefficients of Holder regularity between zero and one arises in a 
variety of geophysical applications. Perhaps the most fundamental one concerns 
the study of thermo-chcmical boundary layers and phase transitions in Earth's 
lowermost mantle — the so-called D" layer overlaying the core-mantle boundary 
(see Figured} cf. [25] for recent images of the phase transformation in D"). Such 
phase transitions can only be probed by earthquake generated seismic waves 
through scattering off these. The relevant scattered wave constituents appear 
as precursors to, for example, the core-reflected compressional PcP phase and 
the horizontally polarized shear ScS phase. The scattering is most prominent 
at large opening angles (towards grazing incidence). It is to this situation that 
the paraxial approximation and coefficient dependence considered here applies. 
In this context the principal direction, z, of propagation is perpendicular to the 
direction of backscattering, x. The central question is whether an imprint of 
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mantle 




Figure 1: The so-called D" layer (in gray) above the core-mantle boundary 
(at approximately 2800km depth) with a core-reflected (ScS) wave and two 
precursors. The precursors in the neighborhood of the top of D" are locally 
modeled by our equation (indicated by ellipsoidal region), x coinciding locally 
with the radial direction and z coinciding with the tangential direction. 

the coefficient regularity on the (regularity of the) scattered wave occurs. A 
positive answer to this question, as provided in this paper, allows to infer a lack 
of lateral regularity in the medium from a lack of regularity of the measured 
data. 

Seismic inverse scattering has been formulated mathematically in terms of evo- 
lution equations with respect to the depth variable in cf. [TTl [HI [19] (with 
smooth symbols in the single scattering approximation) . At the basis of these 
models are one-way wave equations, which are typically of the form 

(d z ±iB(z,x,D t ,D x ))u = f, 

where B is (microlocally) a pseudo differential operator with principal part 

b( Z , X, T, = TC(Z, X)- 1 y/l-T^c(z, X ) 2 \^. 

Here, t G M is time, z > denotes depth, x S M d are lateral spatial directions, 
and D = i~ 1 d. Approximation of the square root to leading order results in the 
standard Schrodinger-type paraxial equation 

1 d 

& + dz~x) Dt ^ W + 2i ^ D *i 4 D *' W ) = °' 

where c(z, x) is the local speed of propagation and z plays the role of the evolu- 
tion parameter. In the frequency (r-)domain the above equation is transformed 
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with a so-called co- moving frame of reference according to w{z, x, r) = u{z, x, r) 
exp(irT(z, x)), where T(z,x) = J Q dz'/c(z',x). Then the paraxial equation at- 
tains the form 

(d z -iA)u = 

in which A is given by 

M = \ E ^™D Xj (c(z, x)r- l D Xi (e— T ^u)). 

The second-order differential operator A can be written as the sum of the self- 
adjoint operator 

1 d 

A(t; z, x, D x ) : = 2 E D *: c ( z ' x)t^ 1 D Xi 
i=i 

and a symmetric perturbation. We observe that if \\D x T(z, •)lli oc (K <i ) ^ ^ 
where K is an appropriate constant, then it is guaranteed that this perturbation 
is ^-bounded with relative bound less than 1 (cf. fTJl Section X.2]). From the 
viewpoint of generators of strongly continuous contraction semigroups on L 2 (R d ) 
the simplification of A by A is of no consequence. 

Remark 1.1 (Directional decomposition and one-way wave equations). 

Directional decomposition leads to the introduction of one-way wave equations 
|21j . One-way wave equations approximate solutions to the wave equation mi- 
crolocally, relative to a principal direction of propagation. (This principal direc- 
tion does not need to be defined globally; one can introduce curvilinear coordi- 
nates and an associated Riemannian metric to generate such directions locally.) 
The validity of one-way wave propagation breaks down when singularities tend 
to propagate in a direction perpendicular to the principal direction (that is, in 
a transverse direction). Hence, to make the statement concerning approxima- 
tion above, precise, one needs to introduce a microlocal attenuation |T6 . The 
mentioned procedures and results require smooth coefficients and symbols, and 
can be proven by making use of the calculus of pseudodifferential operators and 
Fourier integral operators with complex phase. 

It has been demonstrated that, in special cases, one can weaken the condi- 
tion of smooth coefficients. For example, if the coefficients are independent of 
the coordinate along the principal direction, one can allow a step function (in a 
transverse coordinate) and still solve the associated scattering problem by meth- 
ods of one-way wave equations. The approach to carry out such an evaluation 
can be found in [5]. Indeed, scattering in the transverse directions can, at least 
in special cases, be incorporated in the one-way wave equation. Moreover, in the 
case of wave propagation in random media, the (stochastic) paraxial equation 
naturally appears as well [12] . However, a general result concerning directional 
decomposition for - or recomposition to solutions of - the wave equation with 
non-smooth coefficients has not been obtained. 
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In this paper, we address the general problem of "transverse scattering" by a 
one-way wave equation, which constitutes one component in the development 
of a general theory referred to above. 

Remark 1.2 (Regularity and the second order wave equation). Both, 
mode decoupling of a second order wave equation into one-way wave equations as 
well as the derivation of the narrow beam approximation outlined above, require 
higher order differentiablity of the coefficient c(z, x) with respect to all variables 
to make sense (due to truncation of symbol expansions). However, the resulting 
paraxial wave equation displays precisely the same coefficient regularity as in 
the original second order wave equation. Hence as a model equation it still 
reflects the correct medium properties on all scales. The (exact) solutions of the 
paraxial equation then serve as a narrow beam approximation to solutions of the 
original wave equation. In particular, the regularity properties are comparable 
on the same scales. 

The fine tuned well-posedness theorem for wave equations by Colombini and 
Lerner in [3] assumes Log-Lipschitz regularity of the coefficients in the principal 
part. Moreover, their results are sharp in the sense that counterexamples to 
solvability exist when the coefficient regularity is below Log-Lipschitz (but still 
of any continuity type arbitrarily close to such). In order to indicate how their 
results relate to ours, we repeat the key energy estimate from [4]: 

A function a £ L°°(R d ) is said to be a Log-Lipschitz function if 

mi i / \ i \ a ( x ) ~ a (y)\ 

a T T := sup \a[x)\ + sup — — - — < oo. 

LL ie i* \(x - y)log(\x - y\)\ 

|»- v |<l/2 

Colombini and Lerner consider second order wave operators of the form 
Pu := d\u - ^2 (a 3 k(x, t)d Xk u) + M (x, t, d t ,d x )u, 

\<j,k<n 

where (djk)i<j,k<n is re& l symmetric with Log-Lipschitz components and satis- 
fies with some 1 > Sq > the strong ellipticity condition 

]T a jk (x,t)^ k >6o\tf Kel",( Il f)el' 1+1 ) i 

l<j,k<n 

and M(x,t,dt,d x ) is a first-order differential operator with Holder-continuous 
coefficients. 

Let P be as above and 9 s]0, 1/4]. There exist f3 > and C > such that for 
u £ C 00 (M"+ 1 ) and t e [0, 1//3] the energy estimate 

sup \\d t u(.,s)\\ H -e-f3s + sup \\u(.,s)\\ H1 - e -0 S 

0<s<t 0<s<t 

< C (J \\Pu(.,s)\\ H ^- 03 ds+\\d t u(.,0)\\ H -e + IK, 0)11^-*) , 
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holds (cf. [H Equation (2.6)]), where depends only on So, on the Log-Lipschitz 
norm of the djk, and on the Holder norms of the coefficients in M(x, t, dt, d x ). 

We show in the sequel that for the paraxial wave equation the condition on 
the coefficient regularity can be relaxed. For example, if e > any function in 
H 1+£ (M. 2 ) of local behavior like x t— * |2;| 1 / 2+e is not Log-Lipschitz continuous 
but satsifies the assumptions of our main results below. 

The low coefficient regularity in our model conditions has its price in terms 
of a few technical aspects of the current paper: Additional care is needed in 
identifying the appropriate distribution and function spaces that allow for the 
description of mixed regularity properties and for the rigorous formulation of 
a solution concept. The existence proof then consists in showing a series of 
functional analytic properties to establish an evolution system of operators; 
among these the basic self-adjointness property — in the disguise of an elliptic 
regularity lemma — is derived by employing rather delicate regularity properties 
of multiplication in certain subspaces of the space of distributions. 

The plan of the paper is as follows. In Section 2 we present the precise form 
of the operator and specify our (low) regularity assumptions on the coefficients. 
The solution will be sought as a continuous map of depth into the space of tem- 
perate L 2 -valued distributions. Section 3 is devoted to the construction of the 
evolution system in the frequency domain. First, we prove that A(r; z, x, D x ) 
generates a unitary group at fixed z and r. The determination of its domain re- 
quires delicate use of the duality product of distributions as well as a bootstrap 
argument involving multiplication in scales of Sobolev spaces. This leads to 
the construction of an evolution system at fixed r. Finally, strongly continuous 
dependence on the frequency parameter r is established, based on a difference 
approximation. Again a subtle interplay between regularity arguments and dis- 
tributional products is at the heart of our arguments. The strong continuity 
enables us, in Section 4, to construct a solution of the evolution system in fre- 
quency domain with distributional data. Finally, existence, uniqueness, and 
regularity of solutions to the original Cauchy problem is obtained. As an ap- 
plication to inverse regularity analysis, we obtain that a lack of if 2 -regularity 
in the observed solution implies the existence of a region in which the lateral 
regularity of the medium is at most H 1 on the Sobolev scale. 

2 The Cauchy problem: function spaces and 
coefficient regularity 

We recall the definition of temperate distributions on K with values in a Banach 
space E (cf. [531 Section 39.3]; but note that we use a different topology here): 
let J7"(M.; E) be the space of continuous linear maps <y(M) —* E, equipped with 
the topology of pointwise convergence; the Fourier transform T on J?(M) is 
extended to E) by setting (J r G)(0) = G(Tcj>) (G <E J/'QR; E), (f> € >(R)), 
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which is easily seen to be an isomorphism (for the locally convex structure). 

We denote the time variable by t £ M. and introduce coordinates z £ [0, oo) 
for depth (the evolution direction in our context) and x £ R d for the lateral 
variation, where d < 2. As basic space of the wave components we consider 

(2.1) V :=C([0,oc),y'(R;L 2 (R d )). 

Its elements are continuous maps z t— > u(z) of the depth variable z into tem- 
perate distributions of time t valued in L 2 -functions of the lateral variables x. 
When we need to keep track of precise regularity information in the lateral 
variation of the waves, we may employ the Sobolev-scale H s (R d ) (s £ R) and 
define 

(2.2) V s := C([0, oo), S"'(M.; H s {R d )). 

Let T t ■ y'{R; H s {R d )) -» ,Y'(R; H s (R d )) denote (partial) Fourier transform 
with respect to the time variable. We extend Tt to an isomorphism Tt of (the 
locally convex structure of) V s by 

(T t u)(z):=Tt{u{z,.)) Vz>0. 

We consider the following Cauchy problem for a prospective solution u £ V (div 
and grad with respect to x £ R d ): 

(2.3) Pu := d z u - i div (C(z, x, D t ) ■ grad u) = f £ V 

(2.4) u \ z=a - uo £ <y'{R:L 2 (R d )). 

Here, C is a pseudodifferential operator in t with parameters z and x. While in 
the classical paraxial wave equation it is of order 1, here we may assume it is of 
some order m £ R. The precise conditions are collected in the following 

Assumption 1. The symbol of C(z, x, D t ) is of the form 

N 

(2.5) C(z, x, t) = c(z, x, r) • I d = (c a + ^ c l( z > x ) h i( T )) ' h, 

i=i 

where N e N, Id is the d x d identity matrix and the following hold: 

(i) For I = 1, . . . , N : hi is a real-valued smooth symbol (of order m) on R, i.e., 
for all k £Nq an estimate \d^hi(r)\ = 0{\r\ m ~ k ) holds when |r| is large; 
in addition, we assume that \hi{r)\ > rjo near t = (I = 1, . . . ,N) with 
some constant r] > (this can be achieved by adding a cut-off function 
without changing the relevant frequency range). 

(ii) cq is a positive constant. 
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(Hi) there isanr e (0, 1) such that for I < I < N: q is fflC'([0, oo), H r+1 (R d )) 
and real-valued. 

(iv) for all (z, x, r) G [0, oo) x M. d x K; c(z, x, r) > c . 

Remark 2.1. (i) The operator action corresponding to a typical term ci(z, x)hi(D t ) 
in the sum decomposition (|2.5|) on any element w of V is given as follows: for 
all G J^(R) 

which is then multiplied by the function q(z, .). 

(ii) Note that parts (i-iii) of the Assumption imply, for any z and r fixed, the 
following (Zygmund-)Holder-continuity: 

c(z, .,t) -c G F' r+1 (K d ) C Cl +1 ^(R d ) 

(cf. [HI Proposition 8.6.10]). Thus the coefficients have lateral Holder-regularity 
of order r + 1 — |. In the most relevant case from geophysics, d = 2, this yields 
coefficients in C£(R 2 ), but not necessarily Lipschitz continuous. (This includes 
the situation of a boundary layer, as discussed in the introduction, where the 
coefficient is smooth in one of the two variables.) Part (iv) implies uniform 
cllipticity of the lateral differential operator. 

(iii) Note that in dimension d > 3 Sobolev regularity of order r + 1 would not 
imply continuity of the coefficients (in case 1 < r + 1 < d/2). 

Example 2.2. We consider a model with two-dimensional lateral variation 
(d = 2) in the medium properties of low Holder regularity depending on depth. 
In (|2.5[) we put c; = when / > 2 and let c\ be of the form 

Cl (z,x) = X (z,x)\x\ a W 

with the following specifications: x € C 1 ([0, oo) x M. d ) such that x(z, x) = Xo{ z ) 
when \x\ < R\ and x( z > x ) = when \x\ > R2 for certain radii < R\ < R2 
and some xo € C 1 ([0,oo)); a G C 1 ([0,oo)) with some uniform positive lower 
bound oto, i.e., a(z) > a$ > for all z. Then we may choose any r such that 
< r < ao and obtain the following regularity properties at arbitrary fixed 
values of z and r: 

c(z, .,t) - c G H r+1 (R 2 ) n C? {z) (Wl 2 ) C i/ r+1 (M 2 ) n C°°(Wl 2 ). 

(Observe that locally in two dimensions, for any < s < 1, the function \x\ s 
belongs to CI and to H s+1 ~ e for every e > but not to H s+1 .) 

Applying Tt to (|2.3H2.4p we obtain an equivalent formulation of the Cauchy 
problem in the frequency domain: 

(2.6) Pv — d z v — i div (c(z, x, t) grad v) = g G V 

(2.7) v\ z=0 = v Q e^'(R;L 2 (R d )). 
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Equation (|2.6[) is an evolution equation for depth z with the second-order oper- 
ator 

(2.8) A(t\ z, x, D x )v := div (c(z, x, r) gradw) 

acting in the lateral x-domain and smoothly depending on the "external" pa- 
rameter t. Note that A(r; z, x, D x ) is uniformly elliptic by Assumption 1 (iv). 

Remark 2.3. Note that P in (|2.3[) commutes with convolution in the time 
variable. Therefore, damping (or cut-off) of high frequencies in the data of (|2.61 
12. 7[) corresponds to time-smoothing the data in the original problem (|2.3H2.4[) : 
more precisely, if a frequency filter \ G is applied by g :— J-t(f) ■ x( T )? 

uo := ^t(wo) • x( T ) an d w is a solution to (|2.6ti2.7|) then u :— T t (v) solves 

(12.3H2.4p with the data changed to f * x an d «o * X- 

(*) («) 

3 Evolution system 

In this section, we will show that, up to any finite depth Z > 0, the family 
of unbounded operators iA(r;z,x,D x ) (t 6 M, z > 0) generates a strongly 
continuous evolution system (or fundamental solution) {U(t; z\, Z2) ■ < Z\ < 
Z2 < Z} on i 2 (R 2 ) in the sense of [20, Chapter 4] which, in addition, is strongly 
continuous in the frequency variable r 6 R. In a first step we freeze both 
parameters, r as well as z, and construct a strongly continuous (semi-)group of 
operators on L 2 (M. d ). 

3.1 Unitary group at frozen values of r and z 

Notational simplifications: By abuse of notation we will employ the short- 
hand symbols A := A(t; z, x, D x ), c(x) :— c(z,x,t), and c±(x) now denoting 
^ ;>1 q(z, x)hi{r). To summarize, using the above conventions and Assumption 
1, we have 

(3.9) Av — div (c(x) gradv) 

as unbounded, formally self-adjoint operator on L 2 (M. d ) with coefficient 

(3.10) c(x) = c + cx(x) with c > and < ci G H r+1 (R d ). 

We will show that A is a self-adjoint operator with domain D{A) = H 2 (M. d ). 

Remark 3.1. Note that self-adjointness of A with domain H 2 would be im- 
mediate from uniform ellipticity in case the coefficient were smooth. On the 
other hand, self-adjointness on some domain could be obtained in an abstract 
fashion via quadratic forms under mere L°°-assumptions ([HI Section VIII. 6]). 
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However, in accordance with our focus on the interplay of the coefficient regu- 
larity class with qualitative solution properties, we will give an explicit domain 
description in terms of Sobolev spaces, which in addition is uniform with respect 
to t and z. 

Observe that, due to the low coefficient regularity, we also have to establish that 
A is well-defined on all of H 2 (H d ) as an operator into L 2 (R d ). This is included 
in the following lemma as the special case s = 0. 

Lemma 3.2. Let < s < r < 1 and v £ H s+2 (R d ). Then Av £ H s (R d ). 

Proof. Each component of gradu is in H s+1 (R d ). Hence multiplying with 
the i? r+1 -coefficient ci, as well as with the constant Co, is well-defined within 
H s+1 (M. d ) since this is an algebra. □ 

We set D(A) :— iJ 2 (R d ) and note that an integration by parts immediately 
yields that A is symmetric, i.e., D(A) C D(A*) and A* \ D ^ A - ) — A. We proceed 
to show that also D(A*) C D(A) by which self-adjointness will be established. 

By definition, the adjoint operator has domain 

D(A*) = {v e L 2 | for some w £ L 2 : (ip\w) = (Aip\v) for all ij) € H 2 }, 

where ( | ) denotes the inner product in L 2 (R d ). Let v £ D(A*) and choose a 
sequence (ffc)fceN m H 2 (M. d ) which converges to v in L 2 (M. d ). On the one hand, 
there exists w £ L 2 (R d ) such that we have for all test functions ip £ I?(R £i ) 

(A*Vk\ip) = (vk\Aip) — > (v\Aip) = (w\ip) when k — > oo. 

Thus, (A*Vk)k&$ has the distributional limit w £ L 2 (M. d ). On the other hand, 
since c\ £ iJ 1 (R d ) and grad^fc — > gradw in H^ 1 (M. d ) (as k — + oo) we may 
employ the continuous duality product of distributions (cf. [HI Proposition 5.2]) 
H 1 x H^ 1 — > W^l' 1 and deduce that A*v k = Av k — > div (c(x) grad v) = Av 
in Wi~ c ' , hence in the sense of distributions. By uniqueness of distributional 
limits, we deduce that Av — w £ L 2 (M. d ). We obtain 

D(A*) = {v £ 1? | Av £ L 2 }. 

The assertion D(A*) C D(A) = H 2 (R d ) follows now from the following result. 
For later reference, it is stated in slightly more general terms. (We first consider 
the important case d = 2 and leave the case of one-dimensional lateral variation 
for a remark below.) 

Lemma 3.3 (Elliptic regularity). Let < s < r < 1 and v £ H S (M. 2 ) such that 
Av £ H S {R 2 ). Then v belongs to H S+2 {M. 2 ). 

The proof will be based on repeated use of the following three facts, which we 
collect in a preparatory list of "ingredients" : 
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Fact A: Let w 3 G H s j (R 2 ) (j = 1, 2) such that si + s 2 > 0. Then 



w x ■ w 2 S ff So (Il 

where 



so 



min(si, s 2 , si + s 2 — 1) if Si ^ ±1, s 2 7^ ±1, and si + s 2 7^ 

min(si, S2, Si + s 2 — 1 — e) with e > arbitrary, otherwise. 



This is included in the statement of [HI Theorem 8.3.1]. As can be seen 
from the proof therein, one also obtains continuity of the multiplication 
H Sl x H S2 — > H s ° (with respect to the corresponding Sobolev- norms) . 

Fact B: We can find a function F G C°°(R), F(0) = 0, such that 

-^ = - + F( Cl (x)). 
c[x) c 

In particular, we obtain 1/c — l/c = -F'(ci) G i? r+1 (R 2 ). Since r > 0, this 
follows from [21 Theorem 8.5.1] once F is given. To find F, we simply set 
F{y) := —y/(co(co + y)) when y > cq/2 and extend it in a smooth way to 
R such that F(0) = 0. 

Fact C: If v G L 2 (R 2 ) the equation 

(CI) Av — div(cgradu) = A(cv) — div(v grade) 

holds in 2?'(R 2 ), where the occurring products are defined as follows: using 
gradw G H^ 1 we get cgradf G ' by the duality product [TTJ Proposi- 
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tion 5.2]; cv G L 2 since c G L°°; and v grade G L 1 because djC G ff r C L 
Under the stronger assumption v G H r+1 (M. 2 ) we have, in addition, that 

(C2) Av = grad c • grad v + c Av 

in 2?'(R 2 ), with the meaning of the products on the right-hand side as 
follows: since Av G i? r_1 , Fact A applies and yields cAv G H r ~ 1 ; fur- 
thermore, grade and gradv both lie in H r , so that another application of 
Fact A shows that their (Euclidean inner) product belongs to fj nlln ( r > 2r - 1 ) , 

Remark 3.4. Note that formula (C2) represents A as an operator with a 
(Holder-) continuous coefficient in its principal part and Sobolev regularity in 
the lower orders. We observe that in such situation, [8l Theorem 17.1.1] gives 
local solvability in H 2 for right-hand sides in L 2 . However, the latter does not 
imply i? 2 -regularity of any L 2 -solution. For the pure regularity question, it also 
seems that methods based on perturbations of constant coefficient operators do 
not apply either, since A is not necessarily of constant strength (cf. pj Chapter 
XIII]). 
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Proof of Lemma WlA To begin with, we only know that v as well as Av belong 
to H S {R 2 ). The proof proceeds in three steps, successively revealing higher 
regularity. 

Claim 1: v G H s+r (M. 2 ) 

We have grade G H r , so that application of Fact A, noting that r — 1 < 0, 
gives i; grade e H a+r ~ i , hence div(w grade) G H s+r ~ 2 . Since Av G if* we 
deduce from equation (CI) that A(cv) € H s+r ~ 2 , which in turn implies that 
cv e H s+r . Now invoke the decomposition from Fact B and write 

v = — cv + F(c\)cv. 

CO 

The first part clearly is in H s+r , and for the second summand the same is true 
by Fact A. (Note that min(s + r, s + 2i — e) = s + r if < e < r.) 

Claim 2: v G iJ r+1 (R 2 ) 

We may start from v G 7f r+s by claim 1 and proceed inductively to show that 

(*) V G jr+min(l,«+ir/2) j > 0. 

Claim 2 then follows upon choosing j sufficiently large (i.e., j > 2(1 — s)/r steps 
will be required). 

The case j = is just claim 1, so we assume that the assertion holds for some 
j > 0. If s + jr/2 > 1 it is trivially satisfied for larger values of j, therefore we 
assume tj := s + jr/2 < 1 and that 

V G H r+t >. 

Fact A gives V grade G H r+t i ■ H r C ijmin(r,^+2r-l- e ) g i7 min(r,t J +3r/2-l) 

upon choosing e < r/2. Thus, using the short-hand notation rj := min(r + 
l,tj + 3r/2) < r + 1 we may infer that div(wgradc) G H r i~ 2 . Since r j — 2 < 
r - 1 < < s we also have Aw G -ff s C H Tj ~ 2 . Equation (CI) now implies that 
A(cv) G H Tj ~ 2 , hence cv G H Tj . Again by Fact B, combined with Fact A, we 
obtain 

v = —cv + F(ci)cv G H r+1 ■ H r > C H r ' = H r+min ^ +r ' 2 \ 

CO 

which means (*) for j + 1 in place of j. 
CTaim 5: v G K S+2 (R 2 ) 

We use a similar strategy as in the proof of claim 2 and prove inductively that 

(**) V G iJ-m( s +2,l+(i+l)r/2) ^ > L 

Claim 3 then follows when j is chosen sufficiently large (i.e., j > 2(1 + s)/r — 1 
steps are required). 
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The basic case j = 1 corresponds to claim 2, so we proceed with some j > 1, 
under the additional assumption r < Sj := (j + l)r/2 < s + 1 to exclude trivial 
cases. In other words, we assume that 

v e H s i +1 . 

Therefore, again by Fact A and choosing a possibly occurring e < r/2, we 
deduce 

grade • gradt; G H r ■ H s i C #mm(> >Sj +r/2-l) = ff min(r,^) ) 

where we have introduced g 3 := Sj + r/2 — 1. Exploiting equation (C2) and 
noting that Av G i? s , s < r, we extract the information that cAv € H min ( s ' qj K 
Once again we use the decomposition corresponding to Fact B and the statement 
of Fact A to deduce 

Av e H min< * s ' qj) + H r+1 ■ H miu ( s ' qj } C H min< - S ' q ^ 

and a fortiori that Av G H min( - s+2 ^ +2 \ But q 3 = (j + l)r/2 + r/2 - 1 hence 
(/j + 2 = 1 + (j + 2)r/2 and (**) is proven with j + 1 in place of j. □ 

Remark 3.5. The one-dimensional analogue of Lemma l3.3l is more elementary: 
Av = (cv'Y € H S (R) implies cv' € i? s+1 (R) and, since Fact B is valid for d = 1 
as well, we obtain v' G iJ s+1 (R); thus, u G H S+2 (R). 

We summarize the intermediate conclusions from the discussion so far in a 
separate statement, where we appeal to Stone's theorem providing us with the 
exponential unitary group T{z) — exp(izA). 

Theorem 3.6. Let c G C(R d ) satisfy S3.10\) and define Av = div (c(x)gradv) 
with domain D{A) = H 2 (R d ) in L 2 (R d ). Then A is self-adjoint and iA gen- 
erates a strongly continuous unitary group (T(z)) zS k on L 2 (R d ). Moreover, we 
have the following resolvent estimate, valid for A G R \ {0}: 

(3.11) \\(iA - A) _1 w|| L 2 < ^t|t- for all v G i 2 (R d ). 

|A| 

We briefly recall how (T(z)) zg R can be used to construct solutions to the Cauchy 
problem on R d x R 

d z v — iAv = g, v(0) = vo. 
Let v G L 2 (R d ) and g G L 1 (R, L 2 (R d )) then the mi/d solution 

(3.12) v(z):=T(z)v + [ T{z-p)g(p)dp 

Jo 

is in C(R, L 2 (R d )). 

If u G i? 2 (R d ) and g G C(R, H 2 (R d )) or g e C X (R, L 2 (R d )) then w belongs 
to C 1 (R, L 2 (R d )) and is the unique classical solution with pointwise values in 
H 2 (R d ) (cf. [El Section 4.2, Corollaries 2.5 and 2.6]). 
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Remark 3.7. The mild solution (|3.12[) defines a weak solution in the following 
sense: v(0) = v and for all <fi € V(R d+1 ) 

- [ ((v(z)\d z <j>(z,.)-iAcf>(z,.)))dz = [ (g(z)\4>(z,.))dz, 

JR K ' JR 

where ( | ) denotes the inner product in L 2 (R d ). To see this, one approximates 
the mild solution by classical solutions (vk)keN to equations with regularized 
right-hand side and initial data (O Section 4.2, Theorem 2.7]): L 2 -convergence 
of Vk(z) — > v{z) (as k — > oo), uniformly when z varies in compact intervals, 
together with the convergence gu — > g in L 1 (R, L 2 (R d )) implies convergence in 
the integral formula above. 

3.2 Evolution system at fixed frequency r 

Let t € 1 be fixed, but arbitrary. We consider the z-parameterized family of 
unbounded self-adjoint operators in (|2.8[) and put 

A(t;z) := A(t;z,x,D x ) (z > 0). 

Let Z > be arbitrary. We will check that, for every r, (iA(r; z)) z >o defines an 
evolution system (or fundamental solution) {U(r; Zi, 22))z>zi>z 2 >o on L 2 (R d ) 
by applying [20l Section 4.4, Corollary to Theorem 4.4.2, p. 102] (cf. also [13l 
Sections 5.3-5.5]). We have to check that the corresponding hypotheses are 
satisfied. 

First, observe that D(A(t;z)) = H 2 (R d ) is independent of the evolution pa- 
rameter z (and of r), and every iA(r\ z) is the skew-adjoint generator of a 
strongly continuous (unitary) semigroup (T(r, z; C)c>o) on L 2 (R d ). Further- 
more, the resolvent estimates p. lip , valid for all z (and r), immediately imply 
that (A(r; z)) z >q is a stable family of generators with stability constants 1 and 
for all t (cf. M Definition 4.3.1]). 

Finally, we have to check that for all v £ H 2 (B> d ) the map 

[0, oo) 3 z i ► A(t; z)v € L 2 (R d ) 

is continuously differentiable. We may use equation (C2) from Fact C to write 
(with grad taken with respect to x only) 

A(t; z)v = gradc(z, x, r) • grad?; + c(z, x, r) Aw. 

By Assumption l,(ii-iii), we have 

gradc(., ., r) € ^([0, oo), H r (R d )), c(, ., r) - c h Q { T ) G C\[0, oo), H r+1 (R d )). 

Since gradw € H 1 (R d ) and Au G i 2 (M d ) the multiplication rules plus continuity 
properties in Fact A apply (where in case H r ■ H 1 we choose e < r, if d = 2) 
and yield that A(r;.)v G C x ([0, oo), L 2 (R d )). 
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Thus, all hypotheses of [201 Section 4.4, Corollary to Theorem 4.4.2, p. 102] are 
fulfilled. Note that the evolution system is constructed as the strong operator 
limit of discretizations based on the unitary semigroups of each generator, hence 
is contractive. This implies the following intermediate result. 

Proposition 3.8. Let Z > 0. Then for all r € K the family (iA(r; z)) z >q 
defines a unique evolution system (f7(r; Z\, Z2))z>zi>z.2>o on L 2 (M. 2 ) with the 
following properties: The map (21,22) l_ * U(t; Zi, Z2) is strongly continuous, 
U(t;z,z) = I, U(t; z\, z 2 ) is contractive, and 

(3.13) U(t;z 1 ,z 2 )oU(t;z 2 ,z 3 ) = U(t;z 1 ,z 3 ) < z z < z 2 < Z\ < Z; 

moreover, H 2 (M. d ) is invariant under U(t; Zi, z 2 ), for all v S H 2 (M. d ) the map 
(z\,z 2 ) 1— > U(t; Z\, z 2 )v is continuously differentiable, separately in both vari- 
ables, and the following equations hold: 

d 

(3.14) —U(t;z 1} z 2 )v = A(t; Zx)U(t\ z 1 , z 2 )v 

d 

(3.15) - — U(t; zi,z 2 )v = -U(t; z\,z 2 )A{t\ z 2 )v. 



At this stage, we obtain solutions to a version of the Cauchy problem (I2.6H2.7|) 
at fixed frequency r, i.e., 

(3.16) d z v - iA(r; z)v = g e L x ([0, Z], L 2 (R d )) 

(3.17) v\ z=0 = v e L 2 (R d ). 

The mild solution is defined by 

(3.18) v{z) :=U(t; z, 0)v + I U(t; z, p)g(p) dp 

Jo 

and belongs to C([0, Z},L 2 (R d )) ([H Section 5.5, Definition 5.1]). 

Remark 3.9. (i) In the case of classical solutions, we have the following regular- 
ity property: If u G H 2 {R d ) and g G C([0, Z], H 2 (R d )) or g G C 1 ([0, Z],L 2 (R d )) 
then v G C 1 ([0, Z], L 2 (R d )) is the unique H 2 -valued solution and satisfies the 
equation in the strong sense (cf. [131 Section 5.5, Theorems 5.2 and 5.3]). 

(ii) Observe that, at frozen value of r, one may apply [101 Chapter 3, Theorem 
10.1 and Remark 10.2] directly by putting H = L 2 , V = H 1 , and 

d 

a(z; u, v) = y"^(c(z, .)dju\djv) Vu, v G V. 
j'=i 

It suffices to assume c G C 1 ([0, Z], L°°(M. d )), then for any initial value vq G 
H 1 ^) and right-hand side g G L 2 ([0, Z] x R d ) such that d z g G i 2 ([0, Z\,H^ X ) 
there is a unique solution v G C([0, Z],H V ) (1 C 1 ([0,Z],H~' L ) to the Cauchy 
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problem (|3. 16113. l7|) . However, our approach allows for a precise investigation 
of the T-dependence, which is needed to solve the full Cauchy problem (j2.6H2.7p 
with distributional data as well as to transform back to the original problem 
(|2.3H2.4p in Section 4. Furthermore, our results show that lateral 7J 2 -regularity 
of the data is preserved in the solution. 

We thus have established an evolution system in the L 2 -setting. Note that by 
Lemma l3.3l we have, in fact, that A(r; z, x, D x ) is an unbounded operator on H s 
with domain H s+2 for any < s < r. If we were able to establish an evolution 
system on H s then the regularity information encoded into A would be more 
directly preserved. 

3.3 Frequency dependence of the evolution system 

Throughout this subsection, let Z > be arbitrary but fixed. So far, the 
frequency parameter r was arbitrary, but fixed, throughout the construction of 
the evolution system (U(t; Z\, Z2))z>z 1 >z 2 >q- We will prove that the dependence 
on all parameters (r, 2:1,2:2) jointly is strongly continuous. In the sequel, let 
L(E,F) (resp. L(E)) denote the set of bounded linear operators between the 
Banach spaces E and F (resp. on E). 

We begin with an observation on the general level of semigroups and evolution 
systems. 

Lemma 3.10. Assume that 

(3.19) (t, z) ^ A(t; z) is continuous K x [0, 00) -> L(H 2 (R d ), L 2 (R d )) 

(with respect to the operator norm) and 
(3.20) 

(r, z, () 1 ► T(t, 2; £) is strongly continuous R x [0, 00) x [0, 00) — > L(L 2 (R d )), 

where (T(r, z; C))(>o denotes the semi-group generated by A(r; z). Then the map 
(r, 21,22) ► U(t; Zi, Z2) is strongly continuous from B := R x {(21,22): Z > 
21 > z 2 > 0} into L(L 2 (R d )). 

Proof. We inspect the basic construction of the evolution system from the family 
of semigroups in the proof of [13\ Section 5.3, Theorem 3.1] and keep track of 
the additional parameter r € R in our case. For all (t, 21,22) £ B we obtain 
U(t; Z\, 22) as the strong limit of U n (r; z±, 22) (as n — > 00), where U n (r; ., .) is 
the evolution system defined as follows: put z J n = jZ/n (j = 0, . . . , n) then for 
t E R, < y < z < Z let 

U n (r- 2, y) := T(r, z l n ; z - y) if z l n < y < z < z l +\ 
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and 



U n ( T ;z,y):=T( T ,z*;z-z*)> ]J T(r, zi; Z/n) ■ T(r, z l n ; z l n +1 - y) 

l+l<j<k-l 

if z l n <y < 4 +1 < z k n < z < z*+\k > I. 

By (|3.20p the right-hand side of each formula is strongly continuous with respect 
to (t, z, y), and the boundary values, when k = I + 1 and y = z 1 ^ 1 or z = z„ , 
match. Hence U n is strongly continuous on B and ||C/ n (r; z,y)\\ = 1. 

As in [T31 (3.13) on p. 136] we have the following integral representation for the 
action on any v G H 2 

U n (T; z,y)v-U m (T; z,y)v = / f7„(r; z, p)(A n (r; p)- A m {r\ p))U m {r; p,y)v dp, 

Jy 

where A n (r; p) is the piecewise constant approximation of A(t; p) with A„(r; p) := 
A(t;z%), when z\ < p < z* +1 , and A„(t; Z) = A(t;Z). By ([3~T9]) we have 
||^4„(r;p) — A(t; p)\\ L ( H 2 £ 2j — > uniformly for (r, p) in compact sets. Passing 
to the limit m — > oo in the integral representation above yields the estimate 



\\Uu{t; z, y)v - U(t;z, y)v\\ L2 < \\v\\ H2 / || A n (r; p) - A(t; p)\\ L{H 2 , L2) dp. 

By the uniform convergence of A n (r; p) (as n — > oo) we thus obtain (local) 
uniform convergence of J7 n (r; z,y)v, which proves the asserted continuity of 
(r,z,y) ^ U{r;z,y)v. □ 



We have to establish conditions (|3.19H3.20|) in the specific context of the as- 
sumptions described in Section 2. In due course, we will make repeated use of 
(r, z)-parameterized variants of Facts A-C, stated in Subsection 3.1. Note that, 
in particular, the function F used in Fact B does not depend on (t, z). 

Lemma 3.11. If A(t;z) (t € M, z > 0) is defined by H2. 8\) then Assumption 1 
implies Lip schitz- continuity of the map in condition 13.19\) . 

Proof. Let M(t, z) denote the operator of multiplication of pairs (v\, V2) <E H 1 x 
H 1 by the scalar function c(r, z, x) — c(tq, zo, x). We write A(t; z) — A(tq; zq) — 
div o M(t, z) o grad as a composition of operators and get the following norm 
inequality 

\\A(t;z)-A(t ;zo)\\ l{H 2 jL 2 ) 

- Il dlv Wh^xH 1 ,L 2 ) ' W M ( T -' Z )\\L(H 1 y.H 1 ) ' IIS I ' a< ^ \ \ L(H 2 ,H 1 xH 1 ) 

<V2\\M(r,z)\\ L(HlyHl) . 
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To estimate \\M(t, z)(vi,V2)\\ H i xH i it suffices to find an upper bound of 
\\{c{z, .,t) — c(zq, tq)) ^Hjyi f° r v <= H 1 . We have 

(c(z, x, t) - c(z , x, t )) v(x) = 

(l _ro) J grad ( z ' r ) c ( z o + CT ( Z ~ z o), T + ct(t - t )) dcr u(x), 

which, upon taking the -H^-norm with respect to x and assuming max(|z — 
zq\, |t-t |) < 1, yields 

\\(c(z, .,t) - c(z Q , -,to))v\\ h1 < max(|z - z \, |t - t |) 

• sup{\\d z c(z' , .,t')v\\ h1 + \\d T c(z', .,t')\\ h1 ), 

where the supremum is taken over (z',t') € [zq — 1,zq + \] x [to — 1,to + 1]. 
Assumption 1 implies that d z c, d T c both are continuous functions of (z, t) valued 
in H r+1 (R d ), which combined with Fact A gives 

\\d z c^',;T')v\\ Hl + \\d T c(z\.,r')\\ H1 

< C\ \\v\\ m (\\d z c(z\ .,T')\\ m+1 + \\d T c(z', .,r')\\ Hr+1 ) < C 2 \\v\\ H1 

with positive constants C\, C 2 and for all (z', t') € [zq — 1, zo + 1] x [to — 1, to + 1]. 
Combining all estimates we deduce that there is C3 > such that \z — zo \ + |t — 
To I < 1 implies 

\\M(t,z)(v 1 ,v 2 )\\ h i xH i < C 3 max(|z - z \, |t - t |) \\(vi,v 2 )\\ H i yH i, 

which proves the asserted Lipschitz-continuity. □ 

Lemma 3.12. If A(t;z) (t € M, z > 0) is defined by V2. 8\) then Assumption 1 
implies the continuity condition \3.20\) . 

Proof. We apply the Kato- Trotter theorem on convergence of semi-groups (cf. 
[26l Chapter IX, Section 12, Theorem 1]). According to this theorem, we obtain 

T(t, z; C) -> T(t , z ; () strongly as (t, z) —> (t , Zq), 
uniformly on any compact interval containing £, 

thus (I3.20p by uniformity, provided that we show strong continuity of the resol- 
vent map (t, z) (A — iA(r; z))" 1 —: R(X, iA(r; z)) for some A > 0. 

Fix A > and let / <E L 2 (R d ) be arbitrary. Define u(t, z) := R(X,iA(r; z))f e 
H 2 {M. d ), so that u solves 

(3.21) Xu(t,z) -iA(T;z)u(r,z) = f. 

Adding the difference iA(r; z)u(t, z) — i^4(To; zq)u(t, z) yields 

(A - iA(tq; z )) u(t, z) = f + i{A(r; z) - A(t ; zq)) u(t, z) =: f + iw(r, z). 
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Hence, u(t,z) = R(X,iA(T ; z ))(f + iw(r,z)) and it suffices to prove that 
w(t,z) -> in L 2 (R d ) as (t,z) (t ,z ). Applying (C2) from Fact C in 
Subsection 3.1, we may write 

w(t, z) — grad (c(z, x, t) — c(zq, x, to)) • gradu(r, z) 

+ (c(z, x, t) - c(z , x, r )) Au(r, z). 

By Assumption 1, the difference c(z, .,t) — c(zq, .,tq) tends to in i? r+1 (E d ). 
In view of Fact A this implies w(t, z) — > 0, if grad u(t,z) as well as Au(t, z) 
stays bounded. To prove the latter, we take the L-inner product with u on 
both sides of equation (|3.21|) and obtain 

a II u IIl2 - i^{cd X] u\d X] u) = (f\u). 

3 

Note that taking real parts here yields the estimate (|3.11|) , which is ||u|| L2 < 
H/11^2 / A . If we take absolute values of the imaginary parts, we may use the lower 
bound c(z, x,t) > cq and the resolvent estimate to deduce J^j \\d Xj u(T, z)\\ L2 < 
II f\\ L 2 /(Acq), uniformly in (r, z). 

Finally, the boundedness of || Au(r, z)\\ L2 is revealed in several steps. First, note 
that (C2) from Fact C applied to (pOTj) yields 

c(z, ., r) Au(r, z) = i(f — A u(r, z)) — gradc(z, ., r) • gradu(r, z). 

The first term on the right-hand side is bounded in L 2 , uniformly for all (r, z), 
whereas the second term is uniformly bounded in H r ~ 1 by Fact A. Hence 
cAm(t, z) is a bounded family in H r ~ 1 and, combining Facts A and C, we 
find that Att(r, z) is uniformly bounded in H r ~ 1 as well. Therefore, u(t, z) has 
a uniform bound in iJ r+1 -norms. From here we may proceed as in Claim 3 from 
the proof of Proposition [3~3l (with s = 0). Indeed, the arguments used there pre- 
serve uniform boundedness properties throughout, since we have such in H r+1 
already. Thus, u(t, z) is uniformly bounded in H 2 , in particular, Au(r, z) is 
bounded uniformly for all (r, z), which completes the proof. □ 

We summarize the preceding results in the announced continuity statement for 
the evolution system. 

Theorem 3.13. Let (U(t; z\, Z2))z>zi>z 2 >o be the evolution system generated 
by the family of operators A(t; z) (t G M, z > 0), defined in \2. 8\) and sat- 
isfying Assumption 1. Then (r, Z\,Z2) ^ U(r; z\, z%) is strongly continuous 
M x {(zi, z 2 ) : Z > zi > z 2 > 0} L(L 2 (R d )). 

4 Solution of the Cauchy problem 

In this section we present our main results: existence and uniqueness of solutions 
to the Cauchy problem (|2.6H2.7j) in the frequency domain and to (l2.3H2.4p in the 
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time domain. 



If E is a Banach space, let C b (R,E) denote the space of E- valued continuous 
bounded functions. Observe that for any G € C&(R, L 2 ) C ^'(R;L 2 ) the ex- 
pression U(t; zi, z 2 )G(t) is well-defined pointwise for all (r, Zi,z%). Therefore, 
collecting the results obtained so far we arrive at the following assertion. 

Proposition 4.1. If v Q £ C b (R, L 2 (R d )) and g G C b ([0, Z] x M, L 2 (R d )) then the 
formula 

(4.22) v(z, t) := C/(r; z, 0)«o(r) + / U(t; z, p)g(p, r) dp 

Jo 

defines a mild solution v £ ^([0, Z],C b (R, L 2 (R d ))) C V to l[2~oj2~7\ ). Moreover, 
when v is a strong solution then u := v is a strong solution of \2.S\2.J$ with 
initial data uq = J\ vq and right-hand side f = J- i T 1 g. 

For example, the hypotheses leading to strong solvability are satisfied if vq <E 
C b (R,H 2 ) and g £ C b ([0,Z] x R, H 2 ) or g e C x ([0, Z],C&(R, ^ 2 ))- Of course, 
using functions that are bounded and continuous with respect to the frequency 
variable r here is just one simple way to ensure that all constructions described 
above work with all involved objects staying temperate. More generally, it would 
suffice to consider elements in V whose distributional action (with respect to the 
frequency variable) is given by (weak) integration over a continuous function 
(times the test function). 

To apply formula (|4.22p to the original Cauchy problem (|2.3H2.4|) we only need 
to state conditions on the data u$ and / that imply J-fUQ £ C b (M.,H 2 ) and 
Ftf G C b ([0,Z] x R, H 2 ) or T t f £ C x ([0, Z],C b (R, L 2 )). Note that, for example, 
in our physical application such conditions would be met if the source or force 
terms are active only for some finite time interval and vanish otherwise. Then 
by the uniqueness of 7? 2 -valued solutions, as stated in Remark 13.91 we obtain 
the following result. 

Theorem 4.2. Assume that the right-hand side f in equation 12. 3\) satisfies 
either f £C([0,Z},L 1 (R,H 2 (R 2 ))) or f £ C 1 ([0, Z], L X (R, L 2 (R 2 ))) . 

Then for every uq £ L 1 (R, H 2 (R 2 )) the Cauchy problem \2.3\2.1$ has a unique 
strong solution u £ C 1 ([0, Z], =5^'(R, L 2 (R 2 ))) which is H 2 -valued in the following 
sense: for all z £ [0, Z\ and <p £ we have (u(z), <j>) £ H 2 (R 2 ). 

Moreover, u is obtained by inverse partial Fourier transform ( with respect to t) 
of v £ C 1 ([0, Z],C b (R, L 2 (R 2 ))) as defined in Equation \4-22fy , where vo := TtU§, 
g := T t f , and U is the evolution system from Proposition \3.8[ In addition, 
v(z,t) belongs to H 2 (R 2 ) for every (z,t). 

Inverse analysis of medium regularity: We conclude with a brief indica- 
tion of a potential application of Theorem l4.2l to an inverse analysis of medium 
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regularity in wave propagation. Suppose that we are in a model situation where 
parts (i-ii) and (iv) of Assumption 1 in Section 2 are satisfied and the regularity 
property (iii) of the medium is in question; assume that we have the a priori 
knowledge that q S C 1 ([0, oo), L°°(R 2 )) for all I. Let the sources of a seismic 
experiment be calibrated to produce data in accordance with the hypotheses in 
Theorem 14.21 If then the measured wave solution u (or v) fails to display the 
asserted iJ 2 -regularity then we may conclude that for some c\ (and near some 
depth z) there is no r 6 (0, 1) such that the H r+1 -regularity holds; in other 
words, the lateral regularity of the medium there cannot be better than H 1 (on 
the Sobolev scale). If we were in the possession of analogous i? s -results (s > r) 
for the Cauchy problem it would enable us to draw sharper conclusions in such a 
"inverse regularity analysis" . Note that the exact location of the most singular 
region need not be known. For the application hence, precise imaging of the 
singularities is not required prior to the regularity analysis. 
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